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Time  series  variability  in  fluctuating 
ocean  waveguides 

Kevin  D.  LePage 


Executive  Summary:  Sonar  performance  is  adversly  affected  by  the  high 
degree  of  variability  introduced  by  internal  waves  and  other  oceanographic 
phenomena  in  the  ocean.  In  addition,  the  prediction  of  acoustic  propagation  in 
general  is  difficult  because  the  environmental  information  generally  available  in 
databases  is  of  poor  resolution  and  is  highly  uncertain.  The  result  is  that  actual 
sonar  performance  often  varies  significantly  from  predictions.  For  this  reason 
efforts  at  characterizing  the  predictability  of  acoustic  propagation  through  the 
uncertain  ocean  are  required. 

In  this  report  the  variability  of  the  temporal  arrival  structure  of  a  broadband 
acoustic  signal  measured  in  a  fluctuating  ocean  channel  is  derived.  Expressions 
for  the  mean  and  variance  of  the  signal  intensity  are  obtained  as  functions  of 
the  statistical  properties  of  the  sound  speed  uncertainty.  The  ocean  variability 
is  required  to  be  characterized  both  as  a  function  of  depth  and  of  range.  These 
types  of  characterizations  are  often  available  from  gridded  databases  and  are 
also  obtainable  from  CTD  casts.  The  theory  uses  these  characteristics  to  esti- 
mate  the  resulting  statistical  properties  of  the  acoustic  signals.  These  statistical 
properties  can  then  be  used  to  estimate  the  performance  degradation  of  various 
sonar  systems.  As  an  example,  the  performance  degradation  is  estimated  for 
one  of  the  advanced  sonar  concepts  being  developed  here  at  SACLANTCEN. 
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Time  series  variability  in  fluctuating 
ocean  waveguides 

Kevin  D.  LePage 


Abstract:  The  variability  of  signals  propagating  through  an  uncertain 

sound  speed  structure  is  addressed.  Signals  are  assumed  to  travel  in  a  nar¬ 
row  band  adiabatically  in  modes  and  to  experience  fluctuations  in  sound  speed 
which  are  characterized  according  to  the  vertical  and  horizontal  distributions 
of  these  fluctuations.  The  sound  speed  fluctuations  are  assumed  to  affect  only 
the  phase  speed  and  the  group  speed  of  the  modes  in  a  perturbative  way.  The 
changes  in  the  local  phase  and  group  speeds  are  expanded  for  small  perturba¬ 
tions  to  the  sound  speed.  Sound  speed  perturbations  are  described  in  terms  of 
their  statistical  characteristics.  Vertically,  the  sound  speed  fluctuations  are  de¬ 
composed  into  empirical  orthogonal  functions  (EOFs),  while  horizontally  they 
are  assumed  to  be  correlated  on  some  horizontal  length  scale  much  smaller  than 
the  propagation  ranges  of  interest.  Thus  the  cumulative  phase  and  group  speed 
fluctuations  over  the  propagation  path  are  assumed  to  be  distributed  Gaussian 
according  to  the  central  limit  theorem. 

The  framework  outlined  above  is  used  to  derive  the  first  and  second  moments 
of  the  signal  envelope  received  over  an  ensemble  of  ocean  realizations  following 
the  distribution  properties  outlined  above.  Since  the  mean  and  variance  of  the 
expected  signal  are  obtained  in  the  time  domain,  the  stability  of  modal  arrivals 
in  time  can  be  predicted  for  a  variety  of  different  sound  speed  fluctuation  dis¬ 
tributions.  Since  the  phase  and  group  speed  fluctuations  are  linear  in  identical 
terms  involving  the  inner  product  of  the  mode  shape  functions  with  the  EOFs, 
the  fluctuations  of  these  quantities  are  entirely  correlated.  However,  as  the  dif¬ 
ferent  EOFs  express  themselves  differently  on  each  set  of  propagating  modes, 
the  modal  interference  structure  becomes  less  certain  due  to  the  fluctuations. 
The  theory  estimates  this  degree  of  decorrelation  as  a  function  of  the  signal, 
waveguide  and  fluctuation  parameters. 

In  order  to  benchmark  the  theory,  the  first  moment  of  the  short  time  average  of 
the  signal  intensity  is  also  predicted  using  realizations  of  propagation  through 
an  ensemble  of  sound  speed  fluctuations  consistent  with  the  statistical  descrip¬ 
tion.  Excellent  agreement  is  found  between  these  self-consistent  Monte-Carlo 
estimates  of  the  signal  variability  and  the  closed  form  expressions. 

Keywords:  signal  fluctuations  o  empirical  orthogonal  functions  (EOF)  o 

complex  envelope  o  narrowband  o  time  reversal  mirror  o  internal  waves 
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1 

Introduction 


The  fact  that  acoustic  waves  travel  through  significant  fine  scale  and  mesoscale  sound 
speed  structure  between  source  and  receiver  at  basin  and  global  scales  introduces 
significant  variability  between  time  series  which  might  be  predicted  using  only  large 
scale  information  about  the  variation  in  sound  channel  and  actual  received  time 
series.  As  acoustic  waves  travel  through  the  sound  channel,  they  encounter  local 
variations  from  the  mean  sound  speed  structure  predicted  by  data  bases  which  are 
introduced  by  1)  biases  between  database  predictions  of  sound  speed  as  averages 
over  annual  measurements  and  a  particular  realization  at  databases  scales  (typically 
1  degree  resolution),  and  2)  fine  scale  structure  which  is  inadequately  accounted 
for  in  the  databases,  such  as  internal  wave  activity.  Because  of  these  measurement 
uncertainties,  predicted  waveforms  which  coherently  superimpose  modes  across  the 
band  of  interest  to  synthesize  a  highly  coherent  snapshot  of  a  particular  realization  of 
the  sound  field  will  vary  significantly  from  actual  measured  arrival  time  series.  Since 
the  purpose  of  predicting  time  series  is  to  help  understand  what  signals  a  detector 
and  time  of  arrival  estimator  may  reliably  expect,  what  is  required  are  estimates  of 
the  expected  arrival  structure  averaged  over  a  large  number  of  realizations  of  sound 
channel  variability. 
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2 

Theory 


We  propose  a  simple  model  for  the  time  series  received  a  great  distance  from  a  source 
which  is  parameterized  not  only  on  deterministic  variables  such  as  center  frequency, 
bandwidth,  receiver  and  source  depth,  but  also  on  uncertainty  accumulated  along  the 
path  introduced  by  sound  channel  fluctuations.  If  we  assume  that  sound  channel 
perturbations  introduce  local  fluctuations  in  the  wavenumber  and  group  speed  of 
modes,  and  further  assume  that  the  fluctuations  are  not  significant  enough  to  cause 
mode  coupling,  we  arrive  at  the  following  expression  for  the  time  series  received  at 
a  receiver  in  the  limit  of  large  range 


p{t,r,z)  ~ 


^  ^-i(wot-fJk^(r')  dr'-l-fjAkn(r')  dr') 


y/lo  K(r')  dr' 

1  ^ 


(1) 


where  Sn  =  and  A^7i(r'),  A5^(r')  are  the  deviations  in  the  modal  wavenumbers 
and  slownesses  caused  by  perturbations  to  the  local  sound  speed  structure  at  range 
r',  and  the  sub  or  superscripts  o  indicate  that  the  quantity  has  been  evaluated  at 
the  center  frequency.  These  have  been  derived  before  [1,  2]  and  are  reiterated  here 
for  convenience 


AA;n(r') 


Ac(2:,  r') 
'r')p{z,  r') 


dz- 


(2) 


A5„(/)  =  -  (^2  -  Snir'),  ^,  I  dz 

\  k?^ir')J  k^(r')  J-oo  ci{z,r>)p{z,r') 

_9_A  r  Ac(2r,r')  d4>n{z,r')  , 

k°ra{T')  J-oo  cl{z,r')p{z,r')  0a;  "  ’ 


It  is  assumed  that  aU  the  mean  properties  of  the  sound  channel  at  range  r'  have 
been  characterized  in  terms  of  the  modal  properties  Sn  and  (j)n^  and  that  the 
remaining  uncertainties  Afe^  and  A5n  are  well  approximated  by  Eqs.  (2)  and  (3). 
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2 A  Characterization  of  signal  variablity  in  terms  of  Empirical  Orthogonal  Func¬ 
tions 

The  question  of  intermodal  correlation  of  the  deviations  in  the  wavenumber  and 
slowness  to  particular  realizations  in  the  sound  speed  perturbations  Ac  leads  us  to 
adopt  the  decomposition  of  the  sound  speed  fluctuations  into  Empirical  Orthogonal 
Functions  (EOFs,)each  of  which  are  by  definition  uncorrelated  for  the  particular 
data  from  which  they  were  obtained,  but  which,  following  Krolik  [3],  will  also  be  as¬ 
sumed  to  be  statistically  independent.  We  seek  to  understand  how  the  decorrelation 
between  adjacent  modes  dictates  the  expected  behavior  of  a  received  time  series. 
These  decorrelations  are  caused  by  the  differences  in  the  responses  of  the  individ¬ 
ual  modes  to  the  common  sound  speed  perturbations  they  encounter  as  they  travel 
along  the  propagation  path.  The  use  of  EOFs  allows  these  differences  in  response 
to  be  decomposed  into  uncorrelated  components;  if  only  one  EOF  is  used  then  the 
perturbations  between  all  the  modes  remain  entirely  correlated. 

The  sound  speed  perturbations  are  expanded  in  terms  of  range  dependent  coefficients 
on  the  EOF  bases 

E 

Ac{z,r)  =  Y^g^{r)ip^{z), 

e=l 

where  the  bases  are  estimated  from  data  sets  according  to  the  formula 

9{:,e)^V{:,e+l)/Vd^. 

Here  V(:,  e  +  1)  are  the  right  singular  vectors  of  the  singular  value  decomposition  of 
the  sound  speed  profile  covariance 

UtSV  =  ^c(r„)ct(z,r,), 

n 

and  dz  =  D/N  ^  where  D  is  the  depth  of  the  water  column  and  N  is  the  number  of 
receiver  depths. 

The  first  right  singular  vector  V(:,l)  is  proportional  to  the  mean  or  background 
sound  speed  profile  and  thus  is  not  used  to  describe  the  random  component  of  the 
sound  speed  profile.  The  range  dependent  coefficients  are  determined  from  the  EOF 
bases  according  to  the  relation 

g(r)=  (4) 

In  this  way  the  range  integral  of  the  deviations  to  the  modal  wavenumber  and 
slowness  are  determined  in  terms  of  known  constants  and  the  range  integral  of  the 
EOF  coefficients  g 
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where 


and 


rr  ^ 

/  Akr^ir')dr'^Ak^T^en9e 


(5) 


[\Sn{r')  dr'  =  ASnY:^enge 

^0  e^l 


(6) 


Akrj 


4>1{z,t')  dz. 


dr' 


The  range  integrals  of  the  EOF  coefficients  themselves  are  assumed  to  be  zero  mean 
Gaussian  random  variables.  Thus  the  range  integrals  are  assumed  to  cover  many 
correlation  lengths  of  the  individual  ge.  If  the  variance  of  the  individual  is 
then  the  variance  of  the  range  integral  is 

^  lo  lo 

-  (7) 

where  is  the  correlation  scale  of  the  horizontal  variability  of  the  EOF  coefficient 
ge.  If  this  horizontal  variability  is  homogeneous  with  known  power  spectrum 
then  the  variance  may  be  expressed  in  terms  of  the  spectral  integral  [4] 

/CO 

dkPg^{k)sinc^{krl2)^ 

-oo 

where  the  spectral  integral  of  the  power  spectrum  P  with  the  sinc^  function  yields 
the  exact  ratio  of  the  correlation  length  scale  to  the  range  Idr. 
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2.2  The  Short  Time  Average  of  the  received  intensity 


The  Short  Time  Average  (STA)  ^  of  the  square  of  Eq.  (1)  is  given  by  the  expression 


PsTAii>r,z)  ~  27r^^ 


-1=1  m=i  yjjo  k°{r')  dr' k^{r')  dr' 


/AtJ  rAuj 

du>i  / 

-  Afi/  J  — 


Acj 


^  gt  (^^en™AA:em+Ct?i  AS'en"'i^2^‘S'em)Pe 
e—l 


(8) 


where  Aken  =  A5en  =  ASn'^e  ^'iid  Alj  is  the  bandwidth  of  interest. 


Assuming  the  are  uncorrelated,  the  expected  value  of  the  STA  of  the  signal 
intensity  is 


N  N  i  i;  ko{r')-k^ir')  dr' 

{PSTAit,r,z)}  ~  ^  — .  . <f>l{z,)(t>l{z,)(t>l{z)4>l,{z) 


n=i  m=i  Y^/o  K{^')  d-'  k°^{r')  dr' 

/Au/  fAuj 

dui  /  du;2 

~  ^\LtJ  —  /\li) 


n 

e=l 


p-(  AA^en  •“  AA^em +t*^l  A5en -u;2  ASem  r4 /2 


(9) 


Note  in  Eq.  (9)  that  the  deviations  of  the  slowness  and  the  horizontal  wavenumber 
are  entirely  correlated.  This  is  a  consequence  of  the  fact  that  these  deviations  are 
both  derived  under  a  local  approximation  and  are  both  the  result  of  linear  operations 
on  Ac. 


Squaring  the  argument  of  the  exponent  on  the  third  line  of  Eq.  (9)  and  collecting 
like  terms,  we  see  that  the  decorrelation  between  adjacent  modes  at  time  t  and  range 
r  is  accounted  for  through  the  evaluation  of  the  following  two  dimensional  frequency 
integration 


rL\UJ 

/»L,r,tK,Aa;)  =  e"®"-  /  dwj 

J—Au> 

/Act; 

-Au; 


^The  short  time  average  is  taken  over  times  sufficiently  long  to  cancel  out  the  2x  carrier  ripples 
in  the  waveform  which  arise  when  it  is  squared.  Thus  the  short  time  average  is  taken  over  time 
scales  proportional  to  l/f^. 
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where 


0L  =  ^(Afc,\-2A^:enAfcem+AfcL)4</2 

e 

©1  —  ^  “  ^^em)  ^  Sen^ g^'^^e) 

e 

0n  =  J^^Sl^(7ljlel2 

e 

®2  ”  ^  ““  ^^en)  ^ 

e 

0^2  =  X]  ^  Sem(^ljU 

e 

Ql,  =  J2ASL<rlrie/2, 


The  intermodal  decorrelation  term  associated  with  0^^^  is  the  modal  decorrelation 
at  Uo  which  has  been  derived  previously  [3].  The  quadratic  decorrelation  terms  Gn 
and  022  can  lead  to  pulse  spreading  due  to  decorrelation  of  adjacent  frequencies;  the 
way  to  see  this  is  that  in  the  absence  of  the  linear  decorrelation  terms  0i^2)  when 
011,22  a^re  on  the  order  of  one  over  the  bandwidth,  1/  A  a;,  the  frequency  integrals 
are  asymptotically  evaluated  as  the  characteristic  function 

J —Auf 


The  frequency  integrals  in  Eq.  (10)  may  be  performed  in  closed  form  if  instead 
of  modeling  the  signal  as  absolutely  band  limited,  we  adopt  the  use  of  a  Gaussian 
amplitude  spectrum  on  the  signal  of  effective  width  Au;.  Then  Eq.  (10)  may  be 
rewritten  as 


J  —  OO 

J  —  OO 
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Eq.  (11)  may  be  evaluated  in  closed  form.  The  resulting  expression  for  the  STA  of 
the  received  intensity  is  (see  Annex  A  for  the  details  of  the  derivation,) 


{PsTAi^rr,z))  ~  27r^^ 


e^Jo 


-nV'  / 


n=i  m=i  yj Jo  kl{r')  dr' k^{r')  dr' 

^  (®22  +  1/2  A  a»2) 

X  V^/(0ir  + 1/2  A  cu2  -  eyiiei^  +  2/ Au^)) 

X  exp  ( ^2  ~  ~  -  Tm(r))\ 

'^V  4022  +  2/Au;2  J 

y  402i  +  2/Au;2_0y(02^_^l/2Aa;2) 

where  Tn,m{r)  =  Sn,^(r')  dr'. 


rnizs)rMr„{z)rM 


(12) 


Eq.  (12)  evaluates  the  expected  value  of  the  STA  of  the  received  intensity.  It  is  also 
of  interest  to  evaluate  the  fluctuations  of  this  STA  as  a  function  of  sound  channel 
variability.  What  is  desired  is 

<^STA{t,  r,  z)  =  {psTA{i:  z))  -  {PsTa(*>  '■> 
which  may  be  written  explicitly  as 


N  N  N  N 

^“sTA  -  E 

n—1  m-l  n'=l  m'=l 

71  Kir')  dr'J^  k^^,{r')  dr' ko^(r')  dr' f,  k°^.{r')  dr' 


poo 


X 


r  duji  r  dwj  e 

J—00  j—oo 
-00 


/oo  TOO 

duJi^  /  duJ2^ 

-00  J—oo 

E 

X  JJ  exp  {  -{Aken  +  Afcen'  -  (Ak  em  4~  A  S^n  +  ijJ\i  A  5*5 

e=l 


-(^2  A  Sem  +  W2'  A  Sejn')f  /'2^ 
-  {PSTAii,  r,  Z))‘^. 


(13) 


Using  the  same  technique  used  for  the  ensemble  average,  the  second  moment  of  the 
ensemble  average  may  be  evaluated  in  closed  form.  The  result  is  given  in  Annex  B. 
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Narrow  band  approximation  Type  0  vs  Fourier  (1500  Hz,  75  Hz  BW,  10  km) 


Figure  1  Comparison  between  envelope  of  the  Fourier  synthesized  impulse  response 
and  the  envelope  of  the  Type  0  narrow  band  approximation.  Agreement  in  general 
is  poor  because  the  Type  0  narrow  band  approximation  neglects  the  dispersion  of  the 
individual  modal  arrivals. 

Note  that  the  scintillation  index  %  may  be  obtained  from  z)  according  to 

the  relation 

This  measure  is  often  used  to  charcaterize  signal  variability  in  a  random  channel. 


2.3  Dispersive  effects 


In  the  derivations  in  the  previous  section  it  has  been  assumed  that  the  unperturbed 
impulse  response  in  the  waveguide  is  well  approximated  as 


p{t,r,z) 


n=l  VV 

J  —  OO 


(14) 
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Figure  2  Shallow  water  sound  speed  profile  used  in  this  example.  The  profile  is 
downward  refracting  and  overlies  a  slow  isospeed  sediment  5  m  thick. 


which  can  be  evaluated  in  closed  form  as 


p{t,  r,  z)  ~  27r  A  a;  ^ ^ - ^7= - ^s)(i>n{(^o,  z),  (15) 

n=l 


In  reality  Eq.  (15),  which  we  call  the  Type  0  narrow  band  approximation,  is  highly 
inaccurate  since  it  neglects  the  pulse  dispersion  in  the  individual  modes  which  is  so 
important  for  determining  the  coherent  interference  pattern  as  a  function  of  depth 
and  time  at  reasonable  ranges.  This  can  be  seen  in  Figure  1  where  the  envelope 
of  the  Fourier  synthesized  Green’s  function  at  1500  Hz  and  75  Hz  of  bandwidth 
for  a  shallow  water  waveguide  is  compared  to  the  envelope  of  the  Type  0  narrow 
band  approximation.  The  approximation  agrees  poorly  because  it  simply  places 
each  mode  at  an  arrival  time  determined  by  the  modal  slowness  Sn^  neglecting  the 
spreading  of  the  mode  arrivals  in  time  which  leads  to  the  ray-like  interference  effect 
seen  in  the  Fourier  synthesized  result. 


The  sound  speed  profile  for  the  shallow  water  environment  is  shown  in  Figure  2. 
Since  the  slow  sediment  layer  carries  the  first  two  modes  at  1500  Hz,  the  travel 
times  in  this  discussion  are  referred  to  the  first  arriving  waterborne  mode,  Mode  3. 
The  dispersion  curves  for  the  waveguide  are  illustrated  in  Figure  3,  and  show  typical 
shallow  water  behavior.  While  the  phase  speed  of  the  first  two  modes  is  quite  slow 
due  to  the  slow  sediment  layer,  the  phase  speeds  for  modes  3  -*  70  increase  with 
increasing  mode  number,  as  expected.  The  group  speed  of  the  first  two  modes  is 
also  slow  due  to  the  fact  that  they  are  trapped  in  the  slow  sediment,  but  as  expected 
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Figure  3  The  dispersion  curves  for  the  shallow  water  waveguide  at  1500  Hz.  The 
left  panel  is  shows  the  curvature  term  d^kjduj^^  while  the  right  panel  showns  the 
phase  and  group  speed.  70  modes  are  supported^  the  first  two  of  which  live  in  the 
sediment  layer. 


the  waterborne  modes  have  decreasing  group  speed  for  increasing  mode  number.  The 
modal  dispersion  term  is  negative,  indicating  that  at  higher  frequencies  a  given  mode 
is  less  slow  than  it  is  at  lower  frequencies 


Sn{0J  +  Wo)  =  Sn{0Jo)  “  i-Dn|  (w  -  W„)/2, 

which  we  have  also  observed  to  be  typical  shallow  water  behavior.  For  this  reason 
time  series  in  shallow  water  waveguides  arrive  in  an  inverse  frequency  chirp,  with 
the  highest  frequencies  arriving  first. 


2.3.1  Type  1  narrow  band  approximation 


A  better  approximation  to  the  waveguide  Green’s  function  is  obtained  when  the 
second  order  modal  dispersion  term  Dn  =  is  included  in  the  phase  argument 
of  the  complex  exponential 


A"  -i{wot-k^T) 


p{t,r,z)  ~  - 0„(Wo,Zs)?^n(Wo,2) 

I  n=l  V*™’’ 

/OO  ' 

-OO 


(16) 


This  frequency  integral  can  also  be  evaluated  in  closed  form  to  obtain  what  we  call 
a  first  order  narrow  band  approximation 
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Narrow  band  approximation  Type  1  vs  Fourier  {1500  Hz,  75  Hz  BW,  10  km) 


Time  after  Mode  3  arrival  (s) 


Figure  4  The  performance  of  the  Type  1  narrow  band  approximation  is  superior 
to  the  Type  0  shown  in  Figure  1.  The  effects  of  modal  dispersion  are  now  accurately 
modeled,  leading  to  better  estimates  of  the  modal  interference  which  gives  the  true 
coherent  propagation  structure. 


p{t,r,z)  ~ 


3?  (27r)i/2^ 


(17) 


As  shown  in  Figure  4,  the  first  order  narrow  band  approximation  for  the  waveguide 
Green’s  function  does  a  better  job  of  replicating  the  coherent  arrival  structure  of  the 
Fourier  synthesized  result. 


The  importance  of  modal  dispersion  implies  that  a  modal  arrival  time  is  an  abstrac¬ 
tion  which  in  practice  does  not  exist.  Instead,  the  modes  disperse  sufficiently  such 
that  they  always  overlap  their  neighbors  in  time.  This  means  there  is  a  fundamental 
lack  of  resolution  of  the  modal  arrival  time  which  is  introduced  by  the  waveguide 
and  which  in  general  bandwidth  will  not  remove.  The  modal  dispersion  which  al¬ 
lows  the  strong  interference  between  modes  spreads  the  arrival  time  of  an  individual 
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Contributions  of  individual  modes  (20  m  SD,  50  m  RD,  1500  Hz,  75  Hz  BW) 


80 


Time  after  Mode  3  aniva]  (s) 


Figure  5  The  envelopes  of  the  individual  modal  contributions  to  the  total  response 
received  at  a  depth  of  50  m.  This  figure  showns  that  the  estimation  of  modal  arrival 
times  is  not  meaningful  at  this  frequency  and  range.  Individual  mode  envelopes  are 
highly  overlapped  due  to  dispersion  of  the  individual  modes.  The  Type  1  narrow 
band  approximation  shown  in  red  does  a  fairly  good  job  of  predicting  the  response, 
but  under  or  over  estimates  individual  model  arrivals  due  to  the  assumption  that  the 
mode  shapes  are  independent  of  frequency. 


mode  out  over  times  large  with  respect  to  the  time  scales  of  total  summed  response. 
This  effect  is  illustrated  in  Figure  5,  where  the  envelopes  of  the  individual  modal 
arrivals  for  a  receiver  at  a  depth  of  50  m  in  the  water  column  are  shown  for  all  the 
propagating  modes,  along  with  the  envelope  of  the  coherent  sum  which  is  related  to 
what  would  be  measured  in  an  experiment.  The  black  curves  are  the  exact  Fourier 
synthesized  results,  which  extend  up  to  82  modes  over  the  entire  frequency  band  of 
the  synthesis,  while  the  red  curves  are  the  the  Type  1  narrow  band  approximation, 
which  is  obtained  by  assuming  that  the  70  modes  present  at  the  center  frequency 
exist  over  the  entire  frequency  band  of  interest  (instead  there  are  68  modes  at  1425 
Hz  and  75  at  1575  Hz.) 

From  the  results  shown  in  Figure  5  it  is  clear  that  an  experimentally  measured 
quantity  such  as  the  arrival  time  of  a  large  pulse  in  the  total  time  series  has  no 
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2.3,2  Type  2  narrow  band  approximation 


10  20  30  40  50  60  70 

Mode  number 


Figure  7  The  frequency  derivative  of  the  shallow  water  mode  shapes  at  1500  Hz. 
The  mode  shapes  are  rather  strong  functions  of  frequency. 


The  fact  that  the  mode  shapes  can  be  strong  functions  of  frequency  even  in  relatively 
narrow  frequency  bands  compared  to  the  center  frequency  of  interest  implies  that  a 
first  order  Taylor  series  correction  to  the  modes  shapes  of  the  form 


^  p-i{wot-k^r) 


^n(^o  ?  ^5) 


d(l>niz) 


d(jj 


d(l>n  (•2^5) 


£ 


iOo 

,2 


doo 


(18) 


can  be  added  inside  the  brackets  of  Eq.  (16)  to  obtain  improved  performance. 

The  numerical  values  of  dcbj Blj  as  determined  by  finite  dilferences  for  our  example 
waveguide  are  illustrated  as  a  function  of  depth  and  mode  number  in  Figure  7.  Using 
these  values  Eq.  (18)  can  be  evaluated  to  obtain  the  complex  representation  of  the 
correction  which  when  placed  in  the  brackets  of  Eq.  (17)  leads  to  a  second  order  or 
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Narrow  band  approximation  Type  2  vs  Fourier  (1500  Hz,  75  Hz  BW,  10  km) 


Time  after  Mode  3  arrival  (s) 


Figure  8  The  Type  2  narrow  band  approximation  bears  a  strong  resemblance  to 
the  true  channel  impulse  response.  For  this  reason  the  theory  of  time  series  stability 
developed  in  this  report  is  based  on  the  Type  2  narrow  band  approximation. 


Type  2  narrow  band  approximation 

^-i{wot-k^r) 


p{t,r,z)  ~  (2 


(e 


2ir 


V  I  2(Aa;-2  -  iDuV 

i{Snr  -t)  ^  d(t>n{z) 


(t  -  Snvf 


2(  —  iDriv) )  duj 

i{Snr-t)  \  d<f>n{Zs) 


2(Acj“^  —  iDnr) 


(19) 


The  performance  of  the  Type  2  narrow  band  approximation  is  illustrated  in  Figure  8 
superimposed  on  the  Fourier  synthesized  result.  This  agreement  is  considered  to  be 
sufficiently  good  to  faithfully  model  the  overall  characteristics  of  the  time  domain 
Green’s  function  of  a  waveguide.  In  Figure  9  the  constituent  modal  components 
of  the  total  time  series  for  a  receiver  depth  of  50  m  are  illustrated  for  both  the 
Fourier  synthesized  Green’s  function  and  the  Type  2  narrow  band  approximation. 
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Contributions  of  individual  modes  (20  m  SD,  60  m  RD,  1500  Hz,  75  Hz  BW) 


80 


-0.06  0  0.05  0.1  0.15  0.2  0.25 

Time  after  Mode  3  arrival  (s) 

Figure  9  The  performance  of  the  Type  2  narrow  band  approximation  as  a  function 
of  mode  number.  These  results  are  significantly  better  than  the  Type  1  results  in 
Figure  5. 


Comparison  with  Figure  5  shows  the  Type  2  approximation  does  a  better  job  of 
modeling  the  temporal  characteristics  of  the  intermediate  mode  numbers,  allowing 
more  accurate  modeling  of  the  total  response. 


2.4  Inclusion  of  dispersive  effects  in  expressions  for 


The  dispersive  effects  associated  with  the  Type  1  narrow  band  approximation  are 
introduced  into  Eq.  (12)  and  Eq.  (34)  by  making  the  substitution 

1/2  l/2ALo^-iDnr/2. 

This  substitution  neglects  the  effect  of  terms  of  the  form  ADn,  which  are  assumed 
to  be  negligible.  This  approximation  is  generally  found  to  be  valid,  as  will  be  shown 
in  the  results  section.  The  modification  to  Eq.  (12)  to  conform  to  the  Type  2 
narrow  band  approximation  is  more  complicated  and  is  therefore  presented  in  detail 
in  Annex  C. 
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2.5  Time  Reversal  Mirror  performance 


The  theory  described  in  this  section  to  model  waveform  stability  can  be  used  as  the 
basis  for  a  study  on  the  performance  of  a  Time  Reversal  Mirror  (TRM)  [5,  6,  7].  A 
TRM  receives  a  source  pulse  on  a  VLA,  reverses  it  in  time  and  rebroadcasts  it.  At 
and  near  the  spatial  location  of  the  original  source  a  substantial  degree  of  focusing  is 
found  to  occur.  If  we  assume  that  the  TRM  time  series  are  measured  for  an  unper¬ 
turbed  waveguide,  we  can  study  the  effect  of  the  random  sound  speed  fluctuations 
on  the  focus  performance.  Under  the  Type  1  narrow  band  approximation,  the  time 
series  measured  at  a  depth  z  and  range  r  from  a  source  at  zero  range  and  a  depth 
Z5  is 


p{t,z,r\z,) 


4>n{Zs)4>n{z)  f  ^-iu>t+{k+{uj-Wo)S„+(u>-u:of  /2D„)r  -{u-u/o)^ /2Au:^ 

f  ^  ^-iu{t-S„T)^-u,^(l/2Au,^-iDnT/2)  /2q\ 


In  a  TRM,  the  received  time  series  is  reversed  and  re-broadcast  at  the  receiver 
depths.  A  non-causal  time  reversal  operation  is  equivalent  to  a  conjugation  operation 
in  the  frequency  domain 

p{-t)  =  F--^  {rH}. 

Assuming  that  the  re-broadcast  time  series  is  scaled  by  a  factor  of  r/p(z)  to  remove 
the  effects  of  cylindrical  spreading  to  and  from  the  source  position  and  to  enable  the 
invocation  of  modal  orthogonality  in  subsequent  steps,  the  complex  spectrum  of  the 
TRM  sources  is 


p*fuj,z)  =  ^Ill^l}^I!lflg->l>nrg-iwSnr^-u^{l/2AuP+iDnr/2) 


so  that  the  complex  spectrum  received  through  a  perturbed  waveguide  in  the  vicinity 
of  the  source  range  and  depth  is 


p{oj,  Az, 

IV  n  m 


0n( '^5)^771  (-^5  “I" 

y/k!^knp{zN) 


y.^-u;^(-iDm(r-\~Ar)/2)^-u)^(l/2Aiu^-\‘iDnr/2) 


(21) 


In  Eq.  (21)  the  quantities  Akm  and  ASm  indicate  that  sound  speed  perturbations 
have  been  introduced  into  the  waveguide  since  the  measurement  of  the  impulse 
response  p{t,  z^r\zs).  It  is  assumed  that  the  dispersion  Dm  is  unaffected  by  these 
sound  speed  perturbations.  Since  we  are  interested  in  evaluating  the  TRM  refocusing 
performance  in  the  presence  of  water  column  variability  rather  that  performance 
as  a  function  of  the  number  and  distribution  of  the  TRM  sources,  we  make  the 
assumption  that  the  TRM  sources  are  fully  depth  sampling  (even  into  the  bottom 
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for  bottom  penetrating  modes  included  in  the  Type  1  narrow  band  approximation,) 
so  that  the  orthogonality  of  the  modes  may  be  exploited.  We  also  assume  that  the 
mode  shapes  are  not  strong  enough  functions  of  the  water  column  variabihty  for  the 
orthogonality  relation  to  be  seriously  degraded.  Under  these  assumptions  the  sum 
over  N  becomes  an  integral  which  when  evaluated  eliminates  off-diagonal  terms  and 
results  in  the  received  spectrum 

p{iO  Az  At)  =  ^iknAr  ^iAknjr+Ar)  ^iujSnAr  ^iu;ASnir+Ar) 

n 

^^-2^{kn}r ^-^{kn}Ar /2Auj^  Dn^r/2 ^  ^22) 

Discarding  the  unimportant  terms  and  inverse  Fourier 

transform  of  Eq.  (22)  yields  the  complex  time  series  as  a  function  of  time  and 
defocused  range  Ar  and  depth  Az  from  the  original  source  range  and  depth 

_  ~l~  ^iknAr  ^iAkn(T-\-Ar)  ^-2^{kn)T 


Az,  Ar)  =  e 


X  J  (23) 

Since  the  TRM  defined  here  is  non-causal,  t  is  the  defocusing  time  about  zero  time, 
when  the  ideal  TRM  concentrates  its  energy  in  depth  and  range  about  the  original 
source  position. 

The  expected  value  of  the  STA  of  the  square  of  Eq.  (23)  is  analogous  to  Eq.  (9) 

^  ^  et(^n-Ar,„)Arg-2(3?{A:n}+^{A:m})r 

(PTHM(<»Ar,Az))~2x^  - 1,,  lit.  I - 

n=l  m=l 


\kji\  \kj) 
rAu) 


A  /\(lJ 

doji  /  d0J2 
-Au;  J —Auf 

E 


xl[e 

e=l 


"(AA^en  —  Afcem+C^l  ASen  “^2  ASem)^<J^gg  — /2Auf^  /^Au/^ 


(24) 


so  that  a  result  equivalent  to  Eq.  (12)  is  easily  determined  as 

-2(5{fc„}+9{Am}h 


H,  Jt  e‘(fcn-fcm)Ar 


{PTRMit,  Ar,  Az))  ~  2x  XI  Z)  “ijrnr 

n=lm=l 


'>«K.(^s)<P°m(zs)<l>lizs  +  Az)4>^(zs  A  z)y'7r/(0|2  +  l/2  Aw2) 
X  +  1/2  A  _  042/(4012  +  2/ A  a;2)) 

'0|  -  (t  -  Sm  A  r)^  -  i20|(t  -  Sm  Ar) 


X  exp 


X  exp 


40^2 +  2/ A 


0f  -  -  -^n  Ar)]' 

402,  +  2/  A  a;2  -  0^2/(012  +  1/2  A  a;^) 


(25) 
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3 

Results 


3.1  Internal  wave  characteristics  in  shallow  water 


Internal  wave  realizations  were  computed  using  the  internal  wave  realization  gener¬ 
ator  in  the  PROSIM  [8]  simulation  package.  The  internal  wave  realization  generator 
was  executed  using  an  energy  level  of  320  m^cph,  a  default  buoyancy  frequency 
profile  as  given  in  Elliot  et.  al.  [9],  and  a  typical  summer  downward  refracting  shal¬ 
low  water  profile  as  measured  north  of  Formica  in  the  Tyrrhenian  sea  during  the 
PROSIM  ’97  experiment  [10].  The  resulting  perturbations  of  the  sound  speed  from 
the  background  profile  are  illustrated  in  the  left  panel  of  Figure  10.  In  general  the 
results  show  that  the  strongest  perturbations  are  found  in  the  vicinity  of  the  transi¬ 
tion  region  in  the  sound  speed  between  the  iso-velocity  warm  surface  layer  and  the 
colder  deeper  water.  The  first  9  EOF’s  obtained  from  this  realization  are  illustrated 
in  the  right  panel  of  the  same  figure.  The  EOF  mode  shapes  are  in  general  confined 
to  the  upper  third  of  the  water  column,  with  the  shallowest  penetration  for  the 
lowest  order,  and  most  energetic,  modes.  The  EOF  mode  shapes  shown  obey  the 
relation  dz=l. 


£OFmo<i«  $hap»s  for  shallow  wa»rn  GM) 


n*ng»(km)  Mo<te  Shape  ampldide  (1W  nr  +  S  n) 

Figure  10  Internal  wave  realization  and  first  9  EOFs  for  the  shallow  water  envi¬ 
ronment  and  a  spectrum  level  of  1  GM. 

The  standard  deviation  of  the  EOF  power  expressed  in  m/s  is  illustrated  in  the  left 
panel  of  Figure  11.  The  rms  amplitude  of  the  first  EOF  is  approximately  0.6  m/s, 
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Figure  11  RMS  amplitude  and  correlation  length  scale  of  the  shallow  water  EOFs. 


while  for  the  9th  EOF  it  is  about  one  order  of  magnitude  smaller.  The  rms  excursion 
of  an  EOF  as  a  function  of  depth  is  obtained  by  multiplying  the  rms  amplitude  of 
the  mode  function  (units  m/s)  together  with  the  mode  shape  (units  l/x/ST)  scaled 
by  the  square  root  of  the  total  water  column  depth  D  in  meters.  Thus  the  maximum 
rms  excursion  of  the  first  EOF  is  approximately  2.5  jn/s,  while  for  the  9th  EOF  it 
is  only  0.13  m/s.  The  correlation  length  scale  I  of  the  EOFs  is  indicated  in  the 
right  panel  of  Figure  11.  These  results  indicate  that  the  correlation  length  scale  of 
the  EOFs  decreases  with  increasing  mode  number,  with  3500  m  for  the  first 
EOF  and  only  -£9  ~  80  m  for  the  9th.  The  correlation  length  scale  was  estimated  by 
fitting  a  power  law  spectrum  of  the  form 

to  the  spatial  spectrum  of  the  range  dependent  EOF  coefficients 

where  the  range  dependent  coefficients  g  are  determined  according  to  Eq.  (4).  An 
example  of  the  quality  of  one  of  these  fits  is  given  in  Figure  12. 


3.2  Waveform  stability  in  shallow  water 

Using  the  EOF  properties  obtained  from  the  internal  wave  realization  in  the  previous 
section,  the  stability  of  the  arrival  time  envelope  for  the  range  independent  shallow 
water  waveguide  whos  properties  are  illustrated  in  Figures  2  and  3  has  been  com¬ 
puted  using  Eq.  (12).  The  results  are  illustrated  in  Figure  13.  The  top  panel  shows 
the  magnitude  of  the  complex  envelope  of  the  unperturbed  waveguide  response  for 
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rad/km 


Figure  12  The  Direct  Method  power  spectral  estimate  of  the  range  dependent  EOF 
coefficient  of  -tps  and  the  power  law  fit  for  the  correlation  length  scale,  determined  in 
this  case  to  be  138  m. 


the  Type  2  narrow  band  approximation.  The  second  panel  illustrates  the  ensemble 
average  of  the  STA  of  the  intensity  as  estimated  over  50  realizations  of  perturbed 
sound  speed  profiles  using  a  modified  form  of  Eq.  (19), 


{i-Slrf 
2(Au>-2  +  iDlr) 


where 
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Unperturbed  waveform  (1500  Hz,  75  Hz  BW,  10  km) 


Time  after  Mode  3  arrival  (s) 

Figure  13  The  unperturbed  (top),  Monte-Carlo  (middle)  and  closed  form  (bottom) 
time-depth  arrival  structure  for  a  VLA  10  km  away  from  a  source  in  the  shallow 
water  waveguide.  The  uncertainty  introduced  by  the  sound  speed  fluctuations  causes 
a  decrease  in  the  ensemble  average  of  the  intensity,  shown  at  early  time  in  the  lower 
two  panels,  as  compared  to  the  unperturbed  result  in  the  top  panel.  The  agreement 
between  the  Monte-Carlo  and  the  closed  form  result  indicates  that  the  closed  form 
theory  captures  the  variability  without  requiring  Fourier  synthesis  or  Monte-Carlo 
averaging. 

and 

e=l 

The  lowest  panel  is  the  closed  form  prediction  given  by  Eq.  (12).  The  excellent 
agreement  between  the  Monte- Carlo  result  and  the  closed  form  result  indicates  that 
the  omission  of  terms  caused  by  in  the  closed  form  result  have  not  seriously 
affected  its  accuracy  (these  terms  were  included  in  the  Monte-Carlo  result.)  The 
expected  value  of  the  STA  of  the  received  intensity  estimated  by  either  means  as¬ 
sumes  the  absence  of  mode  coupling  and  the  adiabaticity  of  the  sound  propagation 
through  the  sound  speed  perturbations. 

The  results  in  Figure  13  indicate  that  the  sound  speed  perturbations  caused  by 
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internal  wave  activity  have  the  strongest  influence  on  the  lowest  order  modes  which 
contribute  to  the  beginning  of  the  time  series.  This  is  confirmed  by  Figure  14, 
which  illustrates  the  contribution  of  the  SSP  perturbation  EOFs  to  the  various  mode 
numbers  for  unity  EOF  variance  cr^^.  The  results  show  that  the  phase  uncertainty 
Aken  introduced  by  the  EOFs  most  strongly  affects  modes  3  -  40,  with  low  EOF 
indicies  introducing  more  uncertainty  than  higher  ones.  The  slowness  uncertainty 
AS 

en  and  the  dispersion  uncertainty  Den  are  confined  between  modes  15  and  45, 
with  the  higher  EOF  indicies  contributing  more  to  the  uncertainty.  The  fact  that 
the  lower  order  modes  are  strongly  affected  by  the  EOFs  indicates  that  the  earlier 
part  of  the  time  series  is  less  predictable,  as  the  lower  two  panels  in  Figure  13  show. 
Conversely,  since  modal  indicies  between  45  and  70  are  relatively  unaffected  by  any 
of  the  EOFs,  the  latter  part  of  the  time  series  is  shown  to  be  highly  predictable. 


LL 

s 


A  k  for  each  mode  caused  by  the  9  EOFs 


10  20  30  40  50  eo  70 

A  S  for  each  mode  caused  by  the  9  EOFs 


Mode  number 


Figure  14  The  perturbations  Akem  (top),  ASem  (middle)  and  ADem  (bottom), 
caused  by  the  9  EOFs.  The  variability  affects  the  lower  order  modes  much  more 
strongly  than  the  higher  order  modes.  This  explains  the  fact  that  the  ensemble  av¬ 
erage  of  the  intensity  shows  lower  levels  and  more  smoothing  for  the  early  arriving 
low  order  modes,  but  high  predictability  for  the  late  arriving  modes. 
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Unperturbed  TRM  refocusing  at  source  range  (1500  Hz,  150  Hz  BW,  10  km) 


Monte-Carto  TRM  refocusing  (1  GM,  9  EOFs,  1500  Hz,  150  Hz  BW,  10  km) 


-0.02  -0.015  -0.01  -0.005  0  0.005  0.01  0.015  0.02 


Closed  form  TRM  refocusing  (1  GM,  9  EOFs,  1500  Hz,  150  Hz  BW,  10  km) 


-0.02  -0.015  -0.01  -0.005  0  0.005  0.01  0.015  0.02 


Time  (s) 

Figure  15  TRM  defocusing  in  shallow  water.  The  top  plot  shows  the  unperturbed 
focus  performance  in  the  time-depth  plane ^  the  middle  plot  the  Monte- Carlo  average , 
and  the  lower  plot  the  closed  form  solution.  Significant  defocusing  in  the  ensemble 
average  sense  is  expected  in  the  case  where  the  sound  speed  perturbations  are  present. 


3,3  TRM  performance  degradation  in  shadow  water 


In  Figure  15  the  performance  degradation  of  the  shallow  water  TRM  predicted  by 
Eq.  (24)  is  illustrated  in  the  depth-time  plane.  The  top  panel  indicates  the  perfor¬ 
mance  of  the  ideal  TRM,  which  is  weU  defined  in  depth  with  a  temporal  resolution 
St  approximately  equal  to  one  over  the  bandwidth  (in  this  case  ±150  Hz  rather  than 
the  ±75  Hz  used  in  the  previous  examples.)  The  depth  resolution  is  determined  by 
the  number  of  propagating  modes  and  the  unrecovered  loss  of  the  higher  order  modes 
which  the  TRM  cannot  eliminate  (although  a  model  based  TRM  which  artificially 
amplifies  the  higher  order  modes  could  recover  some  of  the  lost  vertical  resolution.) 
In  the  second  panel  of  Figure  15  the  deterioration  of  the  TRM  performance  is  shown 
in  the  Monte-Carlo  estimate  of  the  TRM  depth  time  response  estimated  from  50 
realizations  of  environmental  perturbations  drawn  from  the  distribution  parameters 
derived  in  the  EOF  section.  A  substantial  loss  in  resolution  is  observed,  with  signifi¬ 
cant  depth- time  spreading  along  an  ambiguous  ridge  extending  from  the  true  source 
depth  down  in  space-time  at  an  angle  determined  by  the  particular  physics  of  the 
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Figure  16  Internal  wave  realization  for  the  Munk  profile  (left)  at  1  GM.  The  right 
plot  shows  the  first  6  EOFs  The  perturbations  and  the  mode  shapes  are  confined  to 
the  water  column  above  the  sound  speed  axis. 


propagation  through  the  perturbed  waveguide.  The  lowest  panel  shows  that  the 
loss  in  performance  is  well  modeled  by  the  closed  form  expression  in  Eq.  (24).  The 
results  seem  to  indicate  that  TRM  performance  is  more  strongly  affected  by  waveg¬ 
uide  variability  than  the  time  series  themselves,  which  were  shown  in  the  previous 
section  to  be  less  significantly  affected. 


3.4  Internal  wave  characteristics  in  a  SOFAR  channel 

To  support  the  study  of  waveform  predictability  in  the  deep  ocean,  the  PROSIM 
package  was  used  to  generate  a  1  GM  internal  wave  realization  for  a  deep  water 
Munk  profile  [11].  The  realization  is  illustrated  in  Figure  16  along  with  the  first  6 
EOF  mode  shapes.  The  EOF  mode  shapes  are  confined  to  the  upper  40%  of  the  wa¬ 
ter  column,  with  the  most  energetic  confined  to  the  upper  20%.  The  corresponding 
standard  deviation  in  m/s  and  correlation  length  scales  are  illustrated  in  Figure  17 
The  standard  deviations,  shown  in  the  left  panel,  are  approximately  one  order  of 
magnitude  smaller  than  in  the  shallow  water  case,  while  the  corresponding  correla¬ 
tion  length  scales,  shown  in  the  right  panel,  are  1  —  2  orders  of  magnitude  larger, 
with  the  largest,  for  the  first  EOF,  being  approximately  100  km.  As  explained  in 
the  theory  section,  the  longer  the  correlation  length  scale  £g  the  more  uncertainty  is 
accumulated  along  a  path  (assumption  always  being  that  path  length  is  longer  than 
correlation  length.)  However,  the  uncertainty  accumulated  is  also  proportional  to 
the  power  in  the  EOFs,  which  in  this  case  is  much  lower  than  it  is  in  shallow  water. 

The  EOF  rms  sound  speed  perturbations  for  the  Munk  profile  drop  off  much  faster 
than  for  the  shallow  water  case.  For  this  reason  only  6  EOFs  were  included  in  the 
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analysis.  Even  so  there  is  a  drop  of  almost  one  order  of  magnitude  between  the  rms 
sound  speed  of  the  5th  and  6th  EOF.  The  rms  sound  speed  of  the  6th  EOF  is  only 
0.009  m/s,  while  the  first  EOF  has  an  rms  sound  speed  of  almost  0.13  m/s.  The 
correlation  length  scale  of  the  6th  EOF  is  only  6  km. 


FIMS  sound  speed  vanabiliv  m  SOFAR  EOFs  Coirefation  lenaft  scaJe  ol  sofar  EOFs 


Figure  17  The  RMS  amplitude  and  correlation  length  scale  of  the  EOFs  for  the 
Munk  profile.  The  longest  correlation  length  is  100  km. 

The  evaluation  of  the  internal  wave  realization  for  the  Munk  profile  resulted  in 
sound  speed  perturbations  which  were  noisier  than  the  shallow  water  results.  This 
noisyness  is  not  considered  to  be  endemic  to  the  PROSIM  model,  nor  is  it  considered 
to  be  harmful  to  the  validity  of  the  results  which  are  based  on  this  realization.  The 
effect  of  the  noisyness  is  that  the  EOFs  have  higher  spatial  frequency  oscillations 
superimposed  on  them  which  are  not  considered  to  be  meaningful.  While  in  the 
subsequent  section  these  oscillations  wiU  be  shown  to  introduce  some  noise  into  the 
estimates  of  the  perturbations  of  the  modal  dispersion  for  the  5th  and  6th  EOF, 
the  contributions  of  these  EOF’s  to  the  modal  dispersion  is  so  low,  and  the  overall 
importance  of  the  modal  dispersion  in  general  is  so  negligible,  that  the  results  were 
used  in  their  present  form  without  introducing  errors. 


3.5  Waveform  stability  in  a  SOFAR  channel 


The  Type  2  narrow  band  approximation  was  used  to  estimate  the  impulse  response 
of  the  Munk  (SOFAR)  profile  at  a  receiver  range  of  1000  km  at  a  center  frequency 
of  100  Hz.  The  envelope  of  the  resulting  time  series  for  3  Hz  of  bandwidth  is  shown 
along  with  the  envelope  of  the  Fourier  synthesized  impulse  response  in  Figure  18. 
The  Type  2  narrow  band  approximation  agrees  very  well  for  late  times,  when  the 
axial  mode  arrivals  form  the  SOFAR  crescendo.  At  earlier  times  the  agreement 
quantitatively  is  less  good,  as  the  narrow  band  approximation  fails  to  accurately 
model  the  wavefront  curvature  of  the  arrival  time  series  near  the  turning  point 
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depth.  However,  as  before,  the  agreement  is  considered  to  be  good  enough  that  the 
major  characteristics  of  the  SOFAR  channel  impulse  response  have  been  captured 
by  the  narrow  band  approximation. 


Narrow  band  approximation  (100  Hz,  3  Hz  BW,  1000  km) 


Figure  18  Performance  of  the  Type  2  narrow  band  approximation  for  the  SOFAR 
channel  at  100  Hz  and  3  Hz  of  bandwidth.  The  agreement  is  good  qualitatively  at  early 
times,  and  quantitatively  at  late  times  in  the  proximity  of  the  SOFAR  crescendo. 

The  time- depth  response  of  the  SOFAR  channel  is  in  most  ways  the  opposite  of 
the  shallow  water  channel.  Inspection  of  the  group  speeds  and  modal  dispersion  in 
Figure  19  show  why.  In  the  SOFAR  channel,  the  group  speed  of  the  higher  order 
modes  is  faster  than  the  group  speed  of  the  first  mode,  rather  than  slower  as  in  the 
shallow  water  waveguide.  In  addition,  the  modal  dispersion  is  positive  rather  than 
negative.  The  result  is  that  in  the  SOFAR  waveguide  the  higher  frequencies  are 
slower  than  the  lower  frequencies,  so  that  the  SOFAR  crescendo  is  an  up-chirp.  The 
overall  result  is  that  the  impulse  response  of  the  SOFAR  channel  is  characterized 
by  the  early  arrival  of  high-angle  precursors,  with  a  final  crescendo  on  the  SOFAR 
axis  associated  with  the  lowest  order  modes. 

Using  the  internal  wave  EOF  mode  shapes,  correlation  length  scales  and  power  ob¬ 
tained  in  the  previous  section,  the  ensemble  average  of  the  STA  of  the  received 
intensity  in  a  SOFAR  channel  was  predicted  by  a  Monte-Carlo  average  over  50  re- 
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Figure  19  The  modal  dispersion  d'^k/du’^  (left)  and  the  phase  (red)  and  group 
speed  (blue)  of  the  SOFAR  channel.  Positive  modal  dispersion  implies  that  arrivals 
are  up-chirps  (arrive  low  frequencies  first.) 


alizations  and  through  the  evaluation  of  Eq.  (12).  The  result  is  shown  in  Figure  20. 
The  upper  panel  is  the  unperturbed  impulse  response  from  the  Type  2  narrow  band 
approximation,  the  middle  panel  is  the  Monte-Carlo  average  while  the  lowest  panel 
is  the  closed  form  result.  The  results  show  that  like  the  shadow  water  case,  it  is  the 
lowest  order  modes  which  have  the  most  uncertainty,  reflected  in  lower  levels  for  the 
ensemble  averaged  STA  of  intensity  versus  the  unperturbed  levels  at  late  times  on 
the  SOFAR  axis.  Thus  the  theory  predicts  uncertainty  in  the  crescendo,  but  relative 
stability  of  the  high  angle  precursors.  Inspection  of  the  modal  perturbations  to  the 
wavenumbers,  slownesses  and  modal  dispersion  in  Figure  21  show  the  reason  why. 
As  in  shadow  water,  the  modal  phase,  shown  in  the  upper  panel,  is  most  strongly 
perturbed  by  the  EOFs  in  the  lower  order  modes.  The  slowness  perturbations  and 
the  modal  dispersion  perturbations  are  more  strongly  felt  by  the  higher  order  modes, 
but  the  overad  level  of  these  perturbations  is  apparently  insufficient  to  cause  notice¬ 
able  spreading  in  time-depth  of  the  early  arrivals  On  the  other  hand,  the  phase 
perturbation  is  evidently  sufficiently  strong  that  the  predictability  of  the  crescendo 
structure,  when  many  of  the  lower  order  modes  arrive  at  virtuady  the  same  time,  is 
negatively  affected.  This  is  due  to  the  fact  that  the  presence  at  one  time  of  many 
modes  makes  the  coherent  sum  sensitive  to  the  phase  fluctuations  introduced  by 
the  waveguide  variabidty.  Since  the  fluctuations  are  in  detail  unpredictable,  the 
ensemble  average  has  less  structure  and  lower  peak  values  than  the  unperturbed 
result . 


^The  noise  in  the  estimates  of  /S.Den  perturbation  for  e  —  5,  6,  which  were  determined  by  running 
ORCA  for  the  Munk  profile  with  and  without  an  EOF  mode  shape  perturbation,  is  not  considered 
real.  However,  comparison  between  the  Monte-Carlo  result,  which  includes  this  term,  and  the  closed 
form  result,  which  ignores  it,  showns  that  ADen  in  general  are  inconsequential.  For  this  reason  the 
results  and  conclusions  of  this  section  are  not  affected  by  these  noisy  estimates 
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Unperturbed  waveform  (100  Hz,  3  Hz  BW.  1000  km) 


Figure  20  The  unperturbed  (top),  Monte-Carlo  ensemble  average  of  the  perturbed 
( middle)  and  closed  form  ensemble  average  of  the  perturbed  SOFAR  channel  with 
1  GM  internal  wave  fluctuations.  The  fluctuations  affect  the  predictability  of  the 
SOFAR  crescendo,  reducing  the  expected  value  of  the  intensity  since  the  individual 
modal  arrivals  which  constitute  the  crescendo  are  unpredictable  with  respect  to  one 
another.  The  high  angle  precursors  are  highly  predictable  as  their  propagation  phase 
is  less  strongly  affected. 


While  the  results  shown  in  Figure  20  for  the  predictability  of  a  SOFAR  crescendo 
are  interesting,  they  do  not  reflect  the  true  degree  of  variability  which  is  known  to 
exist  in  long  range  time  series  in  the  real  ocean.  The  spatial,  temporal  and  angu¬ 
lar  characteristics  of  long-range  deep  water  transmissions  have  been  the  subject  of 
significant  study  in  recent  years  [12,  13].  In  general  it  has  been  shown  that  mode 
coupling  effects,  which  are  ignored  in  the  theory  presented  here,  are  important  for 
explaining  most  of  the  variability  in  time  series  observed  from  field  experiments. 
Thus  predictions  for  the  variability  of  the  time-depth  envelope  of  long-range  SO¬ 
FAR  channel  transmissions  obtained  with  the  present  theory  cannot  predict  several 
important  features  commonly  observed  in  experiments,  such  as  the  presence  at  late 
time  of  significant  energy  in  higher  order  modes,  and  the  overall  levels  of  the  variabil¬ 
ity  observed  for  a  fixed  value  of  1  GM  in  an  internal  wave  field  using  a  Munk  sound 
speed  profile.  However,  some  of  the  other  features  observed  in  the  data,  such  as  the 
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Figure  21  The  perturbations  to  the  phase  /\ken  slowness  A^en  (middle) 

and  dispersion  ADen  (bottom)  caused  by  the  internal  wave  sound  speed  perturbation 
EOFs  in  the  SOFAR  channel.  The  phase  uncertainty  is  only  present  for  the  lower 
order  modes.  While  the  higher  order  modes  suffer  some  slowness  uncertainty,  their 
phase  are  sufficiently  insensitive  to  the  sound  speed  fluctuations  so  that  their  coherent 
sum  remains  highly  predictable.  The  noise  in  the  perturbations  to  the  various  modes 
caused  by  the  5  th  and  6  th  EOFS  is  wrong  but  does  not  effect  the  results,  as  explained 
in  the  text. 


relative  stability  of  the  precursor  arrivals,  and  the  relatively  unpredictable  nature  of 
the  final  crescendo,  are  predicted  by  the  theory.  Thus  we  view  the  contribution  of 
the  theory  presented  here  to  be  the  theoretical  corroboration  of  the  relative  stability 
of  the  early  arrivals  under  essentially  a  perturbation  analysis.  However,  the  levels 
of  the  waveform  predictability  estimated  here  are  too  high,  and  do  not  account  for 
the  observations  of  late  arrivals  of  higher  order  modal  energy,  which  only  a  coupled 
mode  treatment  would  predict. 
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Unperturbed  TRM  refocusing  at  source  range  (100  Hz,  3  Hz  BW,  1000  km) 
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Monte-Carlo  TRM  refocusing  (1  GM,  6  EOFs,  100  Hz,  3  Hz  BW,  1000  km) 


Closed  form  TRM  refocusing  (1  GM,  6  EOFs,  100  Hz,  3  Hz  BW,  1000  km) 
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Figure  22  TRM  performance  in  an  unperturbed  SOFAR  channel  (top),  and  the  en¬ 
semble  average  of  the  focus  in  the  time-depth  plane  obtained  by  Monte-Carlo  average 
(middle),  and  as  obtained  by  the  the  closed  form  expression  (bottom).  Performance 
is  more  robust  than  in  the  shallow  water  case,  but  the  environmental  perturbations 
can  still  be  expected  to  open  up  the  depth  focus  from  a  region  less  than  100  m  wide 
to  one  almost  1500  m  wide. 


3.6  TRM  performance  degradation  in  a  SOFAR  channel 

The  degradation  of  TRM  focusing  at  the  1000  km  range  is  illustrated  in  Figure  22. 
The  top  panel  shows  the  refocusing  performance  of  the  unperturbed.  TRM  in  the 
depth-time  plane.  As  in  the  shallow  water  case,  the  resolution  of  the  temporal  focus 
is  dictated  by  the  transmitted  bandwidth  (here  3  Hz),  while  the  depth  resolution 
is  controlled  by  the  number  of  propagating  modes  and  their  loss.  In  the  ORCA 
runs  which  yielded  the  complex  wavenumbers  for  this  study,  only  the  last  20  modes 
suffered  any  propagation  loss  at  all,  due  to  interaction  with  a  lossy  bottom,  and 
only  modes  higher  than  121  suffered  losses  of  more  than  1  dB.  Attenuation  due  to 
relaxation  processes  in  the  water  column  were  ignored.  Thus  the  depth  resolution  is 
quite  good  in  the  unperturbed  case.  The  ensemble  averaged  results  of  the  intensity 
both  for  Monte-Carlo  (middle  panel)  and  closed  form  (bottom  panel)  show  very 
little  decrease  in  the  temporal  performance  of  the  refocusing,  but  there  is  a  sub- 
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stantial  decrease  in  the  depth  performance.  The  focus  region  is  significantly  blurred 
across  the  sound  channel  axis  between  the  depths  of  500  and  2000  meters,  while 
the  unperturbed  focus  region  is  confined  to  depths  between  1250  and  1350  m.  Both 
the  Monte-Carlo  average  and  the  closed  form  solution  agree  on  the  amount  and  the 
spatial  characteristics  of  the  blurring  introduced  by  the  water  column  variability. 
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4 

Conclusions 


A  comprehensive  theory  to  describe  the  stability  of  time  series  received  in  waveguides 
has  been  developed.  The  approach  is  unique  in  that  closed  form  expressions  for 
perturbed  and  unperturbed  impulse  responses  for  waveguides  have  been  derived 
around  a  modal  basis.  Thus  the  theory  offers  a  complimentary  approach  to  ray 
theory  estimates  of  waveguide  impulse  responses. 

The  theory  presented  here  is  useful  for  estimating  the  upper  bound  of  time  se¬ 
ries  stability  or  predictability  for  waveguides  which  experience  random  fluctuations. 
The  theory  requires  these  fluctuations  to  be  characterized  by  some  number  of  un¬ 
correlated  empirical  orthogonal  functions,  each  with  its  own  power  and  horizontal 
correlation  length  scale.  The  theory  estimates  the  degree  to  which  these  sound  speed 
fluctuations  affect  the  phase  of  the  individual  modes,  and  the  time  of  flight  of  the 
mode  packets  to  the  receiver.  The  more  EOFs  are  used  to  describe  the  water  column 
variability,  the  more  the  acoustic  modes  become  uncorrelated  from  one  another.  As 
the  modes  become  uncorrelated,  their  coherent  sum  becomes  less  predictable,  and 
consequently  the  ensemble  average  of  the  time  series  becomes  smoother  and  lower 
in  amplitude  than  any  particular  perturbed  or  unperturbed  realization.  Thus  the 
theory  predicts  the  “defocusing”  of  the  arrival  structure  in  the  time-depth  plane 
which  occurs  over  many  independent  realizations.  In  practice,  this  defocusing  can 
be  expected  to  occur  over  time  scales  of  the  ocean  process  which  introduces  the 
sound  speed  fluctuations. 

The  theory  estimates  the  upper  bound  of  the  time  series  predictability  because  it 
relies  on  the  least  harsh  model  for  propagation  through  water  column  variability; 
the  adiabatic  approximation.  While  the  mode  shapes  and  the  wavenumbers  are 
acknowledged  to  be  strong  functions  of  frequency  in  the  underlying  propagation 
theory,  the  sound  speed  perturbations  are  assumed  to  affect  only  the  travel  time  and 
the  accumulated  phase  of  the  individual  modes.  Thus  mode  coupling,  the  exchange 
of  energy  into  and  out  of  modes  caused  by  abrupt  transitions  of  sound  speed  or 
even  by  long  accumulations  of  small  coupling  effects,  is  ignored.  Even  more  prosaic 
effects  such  the  change  in  the  mode  shapes  due  to  the  sound  speed  perturbations 
are  neglected.  These  assumptions  imply  that  the  theory  takes  a  view  of  the  effects 
of  water  column  variability  on  modal  propagation  which  is  almost  benign. 

The  simplifying  assumptions  make  possible  the  derivation  of  closed  form  expressions 


-33- 


SaCLANTCEN  SR-319 


for  the  ensemble  average  of  the  STA  of  the  intensity  received  in  a  waveguide.  For 
the  first  time,  these  types  of  expressions  are  now  available  for  evaluating  time  series 
stabiUty  as  a  function  of  the  modal  characteristics  of  the  waveguide  at  the  center 
frequency  and  the  spatial  characteristics  of  the  sound  speed  perturbations  in  the 
waveguide.  The  requirements  for  running  Monte-Carlo  simulations  are  removed, 
and  Fourier  synthesis  is  no  longer  required,  so  that  a  lower  bound  analysis  of  the 
magnitude  of  the  variability  effects  can  be  quickly  obtained  in  closed  form  with  the 
insight  that  affords. 

The  theory  was  used  to  estimate  the  effects  of  internal  wave  activity  on  the  pre¬ 
dictability  of  time  series  in  shallow  water  and  in  a  SOFAR  channel.  The  results 
showed  that  for  reabstic  sound  speed  perturbations  caused  by  a  Garrett-Munk  spec¬ 
trum  with  an  amplitude  of  320  m^cph,  the  lowest  order  modes  in  both  propagation 
scenarios  had  the  most  unpredictable  phases  and  arrival  times.  In  the  shallow  water 
profile  the  result  was  that  the  onset  of  the  arrival  structure,  associated  with  the  low¬ 
est  order  modes,  became  less  predictable,  while  the  later  arrivals  associated  with  the 
higher  order  modes  and  bearing  a  strong  relation  to  bottom-surface  multiples  were 
quite  stable.  In  the  SOFAR  channel  the  crescendo  on  the  channel  axis,  again  asso¬ 
ciated  with  the  lower  order  modes,  was  the  least  predictable  component  of  the  time 
series,  while  the  higher  angle  precursors,  while  somewhat  unpredictable  in  travel 
time,  none  the  less  proved  themselves  to  be  relatively  insensitive  to  the  sound  speed 
perturbations  introduced  by  1  GM. 

Finally,  the  theory  was  utilized  in  its  first  application,  to  estimate  in  closed  form 
the  expected  value  of  the  defocusing  in  the  ensemble  sense  of  a  time  reversal  mirror. 
The  results  shown  for  the  shallow  water  case  indicated  that  1  GM  could  introduce 
substantial  defocusing  in  time  and  depth  even  while  the  time  series  themselves  ap¬ 
peared  to  be  quite  predictable  and  stable.  For  SOFAR  propagation,  the  performance 
degradation  was  less  severe,  but  both  results  indicated  that  once  the  correlation  time 
scale  of  the  water  column  is  exceeded,  significant  blurring  about  the  focus  depth, 
and  by  extension  focus  range,  can  be  expected  to  occur.  While  this  was  known 
before,  a  lower  bound  on  the  defocusing  may  now  be  quickly  and  easily  obtained. 
Since  the  results  are  given  in  closed  form,  it  is  hoped  that  a  better  understanding 
of  the  defocusing  phenomena  may  be  formed. 
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Annex  A 

Evaluation  of  frequency  integrals  for  (psta) 


We  wish  to  evaluate  the  following  frequency  integrals 


P nm ,r o ^ 


J—oo 

/oo 

di^2  eXp(-U;2(02  -  iTm  -  ‘<^1012)) 

-OO  '  v  '" 


B2 


exp(-u;J(0^2  +  1/2  A  a;^)), 


(27) 


A2 


where  =  ^  —  fo  dr'  for  compactness  of  notation.  Completing  the  square, 

the  U2  integral  may  be  evaluated,  reducing  Eq.  (27)  to 


/OO 

dui  exp(-a;i(0f  +  ir„  -  (210^27-,^  -  20^0f2)/4A2)) 

-OO  ^ - V - ^ 

•  exp(-u;i  (0fi  +  1/2  A  o;^  -  012/4A2)).  (28) 

V - ' 


Completion  of  the  square  to  the  argument  of  the  exponential  inside  the  integral  over 
0^1  and  evaluation  of  the  integral  yields 


/>L,r,tK,Aw)  =  sjrjAxyj-KlA^e  ®"”> 

/4Ai  ^(02 -iTmYlAA2  ^ 


(29) 


which  upon  explicit  evaluation  of  the  parameters  Ai,  Bi  and  A2  yields  Eq.  (12)  as 
desired. 


SaCLANTCEN  SR-319 


Annex  B 

Evaluation  of  frequency  integrals  for  (t^ta 


Upon  squaring  out  the  argument  to  the  probability  term  on  the  fifth  line  of  Eq.  (13), 
we  find  it  useful  to  define  the  following  constants  for  the  evaluation  of  the  frequency 
integrals  over  cji,  s^nd  oj2f 


nmn*m' 

AAjgTT,  "h  A/ugn'  A/jg772,  AAjgy72')<^^g^^e/2 

e 

el 

= 

^  ^(AA*en  +  AAjgyi/  —  AA^gTTi  —  ^kem')  A  Sen^g^'^^e) 

02 

—  ^  ^(AfegTT,  -f"  AAlgri'  “  Afcgyyi  —  A  Sem^g^'^^e 

0?' 

^"^(Afcgrt  A  AfcgT^/  —  AfcgTTT,  ^  /S,kQjjit)  A  5'en'^^e^^e) 

el 

—  ^  ^(Afegyt  +  /\ken'  ^k^TTi  —  AA^eTTi^)  A  S Qrn*(7 

e^ 

5;]A5,V,VV2  +  l/2Aa;2 

012 

= 

^  A5e„  A  Sem<^lj^e 

e 

0?i' 

= 

A5en  A  5en'C^p 

e 

0?2' 

]^A5e„  A5'e,„-(T|^r4 

e 

012 

= 

^A5L4r4/2+l/2Au;2 

0ii' 

^  ^  ASgTTi  A  Sen'^g^'^^e 

e 

022' 

^  ^  ^Sern  A  S^rn*^ g^'^^e 
e 

0?'i' 

^  A5e„-(r|^r4/2  +  1/2  A 

e 

0^2' 

^  ^  ASgyj,/  A  S^rjifcrg^vt^ 

e 

02'2' 

X)  A5eV<^5e''4/2  +  1/2  A 

e 


Using  these  constants,  the  frequency  integrals  may  be  rewritten  in  the  condensed 
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format 

Pnmn'm'  =  ^ 


e“®L„'m'  f°°  dull 
J-oc 

/oo 

“OO 

/oo 

duv  e-‘^I'(®^-‘^-’®^2+‘^l®^l+’V('•))e-‘^^©^l, 

•OO 

/OO 

duj2’  exp(-a;2'  {Ql,  -  u;i/02»i/  +  a;202/2  -  “  »7'm'(»'))) 

-OO  l■■..llll■v>  .1.  I  ✓ 


scp(-u;|/  02/2')^ 


where  henceforth  Tn^m  =  Integration  over  L02f  yields  a  term  of  the  form 


yj7r/A2^  exp{Bl/4Ar). 

Squaring  out  82'  and  separating  out  terms  which  are  independent  of  u;i,  uj2  and  Ui*, 
Eq.  (30)  may  be  written 

Pnmn^m*  =  yTT/M^  eXp(C7|//4 

/oo 

duji  exp(-u;i(02  +  zrn(r)  +  0^q2(0^,  --  i>mO/402^2O) 

-00  ^  yi  I 

^12' 

exp(—a;j(0jj  — 02/1/402/2')) 

- - V - ' 

/oo 

dw2  exp(-u;2(02  -  l>m(r)  -20^,2(02/  -  l>m')/40L')) 

-oo  .  ✓ 

C22/ 

•  exp(a;2(wi(0^2  “  2©2'202'i/40L'))) 

C21 

exp(— 0^2  (  022  “®2'2/'^®2^20) 


X  /  du>i>  exp(-wi-(0?,  +  ir„-(r)  +  20|,i,(0^,  - 

J  —  OO  — .  . 

Ci- 

exp(a;i,(w2(0?,2  -  20^,1,02,2/40^,2'))) 

N— - - - ^ 

C21/ 

exp(-a;i-(wi(0?,i  -  20^,i,0^,i/40^,2,))) 

- V - ' 

Cn’ 

exp(— a;i,(0i,i,  —  02'i'/4/©2'2')),  ( 
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(^2'  —  ' 

Completing  the  square  of  the  argument  to  the  exponential  kernel  in  the  wj/  integral, 
we  obtain 

Pnmn'm'  =  y^Tr/Any^Tr/Aj' exp(Ci2,/4AiO  exp(C'|,/4^2')e"®""''-'’"' 

/OO 

du^i  exp(-mi(0i  +  «V„(r)  +  <^12'  -  2(7ii<Ciy4^i')) 

-OO 

•  exp(-wf(0fi  +  Ci2'  -  C111/4A1')) 


,  02  -  j7Vn(r)  +  C22' +  2C'21'C'i'/4/1i' 

xJ_J^2exp  -U2  +^i(2C'2i-C'n'/4Ai.-C'2i) 


cp(-~a;|(022  +  C2'2  ~ 


If  we  define 


(^2  =  ©2  —  iT-m  +  C22'  +  2C2vCif  /4:Ai/ 


C12  =  2C2vCiif/4:Aii  —  C21, 

then  the  U2  integral  yields 

^  7r/A2  exp((C|  +  u;i2C2Ci2  +  C'12)/ 4^2)* 

Then  Eq.  (32)  may  be  written 


^jK/A2^ir/Ary/Tr/A2<exp{Cl/4:A2)exp{Cl/4Ai,)exp{Cl/4A2')e  ®Wm' 

X  f  dui  exp(— u;i(0i  +  i'rn(r)  +  (7i2'  —  2Cii/Ci//4Ai/  —  2Ci2C'2/4A2)) 
y  _oo  ^  ^ 

Bi 

exp(-a;?(02,  +  “  CIJ4A2)),  (3J 


Evaluation  of  the  uji  integral  yields  the  desired  solution 


pnmn'm'  = 

X  exp(C^l4Ai)  exp(C|/4A2)  exp(C'f,/4Ai/)  exp(C|,/4A2'),  (34) 


where 


Cl  =  B^/4Ai. 
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Annex  C 

Evaluation  of  frequency  integrals  for  Type  2 
correction  to  expressions  for  {psta) 


For  PsTA  Type  2  narrow  band  approximation,  we  wish  to  evaluate  the  follow 
integrals  (assuming  for  the  moment  that  a  single  mode  function  is  a  function  of 
frequency,  this  simplifies  the  expressions  but  does  not  affect  the  final  answer  as 
the  necessary  cross  terms  in  terms  of  the  expressions  derived  here  can  easily  be 
determined”) 

J  diJi  +  Wi^) 

X  J  <kj2  (4>m  +  W2 ^^)  exp ( -UJ2  (©2  -  tVm  -  ^1012))  exp( (022  +  T^)) ■ 

B2  A2 

But 

so  what  remains  to  be  evaluated  is  the  integral 

j  (Lit 

X  y^eS|<".)/«.  ,  (35) 

which  may  be  written 

J  (kji  (<^n  +wi^) 

X  (36) 

which  may  be  written  as  four  integrals,  the  first  of  which  is  a  normal  Gaussian 
integral,  the  second  and  third  have  the  form  of  a  first  moment,  and  the  third  the 
form  of  a  second  moment 

dwi  (<^„  (ym  - 
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+  /i,.( 


0ii-2’')e®2/4A2 


2  -0,1(02 +irn)  -a,?(02j-T2)gB|/4A2  I  (37^ 

duj  du  2A2  ^  f  •  V  ^ 


Eq.  (38''  may  be  rewritten 


y^{/  ‘^‘"1  -  2^1:7 '  ’ " 


+  /  dui 


/ 

/ 


2A2 


_  ^-o,iBig-a,2Ai^B|/4A2 

\ 

Wl  e-‘^iSig-o,2AigB|/4A2 


2A2  du  J  du 


+  / 


which  equals 


\ 


d<j>n  d(j>m  Qir 

d(jj  du)  2A2 

\  i  '} 


I 


.2  -a,iBig-u,jAi  B|/4A2 


I2A2  ^  (38) 


IA2  /Ai  e(®2  -"'")"/4A2 

where  and  Ai  are  defined  in  the  Annex  A. 

We  know  that 

y  („.  +  |l)  ,  0, 

SO  the  second  line  of  Eq.  (39)  equals 


e(®i-‘^-)V4A2  gBf/4Ai 


(39) 


For  the  second  moment  we  know  that 
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SO  the  third  line  of  Eq.  (^39)  equals 


Thus  the  expression  for  has  as  its  kernel  the  expression 


X 


4^71  (  4^m 


B'2 


d<t>n 


2A2  du> 


(.  d4>„.  Qh  ,  (. 

du  2A2  ^ 


2A2  du  ) 


/J_  , 

\du}  du  2A2)  \2Ai  ^4^1// 


d(l>n\  Bi 

du)  I  2Ai 


0- 


> 


# 
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